Using Off-diagonal Confinement as a Cooling Method 
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In a recent letter [Phys. Rev. Lett. 104, 167201 (2010)] we proposed a new confining method for 
ultracold atoms on optical lattices, which is based on off-diagonal confinement (ODC). This method 
was shown to have distinct advantages over the conventional diagonal confinement (DC), that makes 
use of a trapping potential, such as the existence of pure Mott phases and highly populated conden- 
sates. In this manuscript we show that the ODC method can also lead to lower temperatures than 
the DC method for a wide range of control parameters. Using exact diagonalization we determine 
this range of parameters for the hard-core case; then we extend our results to the soft-core case 
by performing quantum Monte Carlo (QMC) simulations for both DC and ODC systems at fixed 
temperature, and analyzing the corresponding entropies. We also propose a method for measuring 
the entropy in QMC simulations. 
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I. INTRODUCTION 

With recent experimental developments on cold atoms 
in optical lattices, the interest in the bosonic Hubbard 
model P, IH has dramatically increased. This model 
is characterized by a superfluid-to-Mott quantum phase 
transition for large onsite repulsion and integer values 
of the density of particles. In actual experiments the 
atoms are confined to prevent them from leaking out 
of the lattice. This is currently achieved by applying 
a spatially dependent magnetic field Q . A parabolic po- 
tential is added into the Hubbard model [|[ to mimic 
the effect of the magnetic field. Therefore, the resulting 
model does not exhibit a true superfluid-to-Mott tran- 
sition, since Mott regions always coexist with superfluid 
regions. This was predicted theoretically Q, and later 
confirmed experimentally [f|. 

Recently, we have proposed a new confining tech- 
nique Q where the atoms are confined via a hopping 
integral that decreases as a function of the distance from 
the center of the lattice. Since the confinement of the par- 
ticles is due to the hopping or off-diagonal operators, we 
called it Off-Diagonal Confinement (ODC), as opposed 
to the conventional diagonal confinement (DC) which 
makes use of a parabolic confinement potential that is 
reflected in the density profile 0| • For large on-site repul- 
sion the ODC model exhibits pure Mott phases at com- 
mensurate filling while at other fillings it exhibits more 
populated condensates than the DC model. Another ad- 
vantage of ODC is that simple energy measurements can 
provide insights into the Mott gap, while the presence of 
the harmonic potential may renormalize the value of the 
gap with respect to the uniform case Q- 

In this paper, we show that the ODC method can also 
lead to lower temperatures than the DC method for a 
wide range of parameters. Producing low temperatures 
in experiments is challenging, especially with fcrmions for 
which laser cooling is not as efficient as for bosons. In 
current experiments, fermions are cooled down by con- 
vection in the presence of cold bosons, leading to Bose- 



Fermi mixtures [§4TlT| . Achieving lower temperatures for 
bosonic condensates will therefore result in colder Bose- 
Fermi mixtures. 

This manuscript is organized as follows. In section II 
we define our model and describe our methods. The hard- 
core limit is studied in section III in order to illustrate 
analytically how ODC produces temperatures that are 
lower than those obtained with DC. This will also serve 
to benchmark the quantum Monte Carlo (QMC) simu- 
lations we use for analyzing the general soft-core case. 
In section IV we present the algorithm we use for QMC 
simulations, and we propose a method for measuring the 
entropy with this algorithm. Results for the soft-core 
case are presented in section V. Finally we conclude in 
section VI. 



II. MODEL AND METHOD 

We consider bosons confined to a one-dimensional op- 
tical lattice with L sites and lattice constant a = 1. The 
Hamiltonian takes the form: 
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The creation and annihilation operators a\ and a i sat- 
isfy bosonic commutation rules, [o,-jOj] = [ a ]\ a j] = 0, 
[a.^aj] = Sij, and hi = a\a i is the number of bosons 
on site i. The sum 53 (jj) runs over all distinct pairs of 
first neighboring sites and is the hopping integral 
between i and j. The parameter U is the strength of the 
local on-site interaction, and W describes the curvature 
of the external trapping potential. 

In this work we consider the grand-canonical partition 
function, 



Z = Tr e 
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where j3 
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, ks is the Boltzmann constant and T 



the temperature. The chemical potential /i controls the 
average number of particles, N = (JV), with N = J^i 
The conventional DC model is obtained by setting Uj = 1 
for all pairs of first neighboring sites and using W > 
0. For this model the value of L is irrelevant as long as it 
is sufficiently large to contain the whole gas. The ODC 
model is obtained by setting W — and using a hopping 
integral Uj that decreases as a function of the distance 
from the center of the lattice, and vanishes at the edges. 
For this model, L fully determines Uj as described below. 

Typically the temperature is not a control parameter in 
cold-atoms experiments, and once laser cooling has been 
performed, the system has a fixed entropy which can be 
considered the control parameter. Then the temperature 
can be estimated numerically knowing the isentropies of 
the system 12]. Therefore, our strategy for determining 
which of the two confining methods can achieve the low- 
est temperature is based on switching adiabatically from 
DC to ODC, so the entropy is conserved. Then we deter- 
mine the temperatures Tdc and T D d c of the DC and ODC 
systems by equating the entropies 

We will consider an experiment in which a fixed num- 
ber iV of atoms is loaded into an optical lattice with a 
DC trap, described by Eq. (JlJ with parameters Uj = 1, 
W = 0.008 (as in Ref. [|). We use L = 400 in order 
to ensure the confinement of the whole gas. Then, we 
adiabatically switch to the ODC trap by slowly varying 
Uj and W to t tJ = (i + j + 1)(2L - i - j - 1)/L 2 with 
L = 70, and W = (as in Ref. 0), keeping N and U the 
same. 

However, in our calculation, it is actually more conve- 
nient to control the temperature than the entropy. Thus 
we consider both DC and ODC systems for a set tempera- 
tures T, and measure the corresponding entropies, Sd c (T) 
and S dc(T). Then, knowing the initial temperature Td c , 
the final temperature T a d c can be extracted graphically 
by imposing the equality, S dc {Tdc) = S dc(T dc), as de- 
scribed in the next section. 



III. THE HARD-CORE CASE: EXACT 
ANALYTICAL RESULTS 

The hard-core limit (U = +oo) of the model can be 
solved analytically. These exact results provide a solid 
benchmark for our study of the general soft-core case in 
the next section. We follow here the method used by 
Rigol [13]. In the hard-core limit, the U term in (JTJ) can 
be dropped if the standard bosonic commutation rules 
are replaced by [a^a-] = [ a L a ]] = [ a i' a j] = f° r 
i ^ j, and a 4 a 



- a] a. 



1, and a,- 2 = a] 2 — 0. With this 



algebra, the model ([T]) reduces to 



# = - E * 



, n alaj + h.c. 



L/2) 2 h u (3) 



which describes hard-core bosons. By performing a 
Jordan- Wigner transformation, the hard-core creation 
and annihilation operators can be mapped onto fermionic 
creation and annihilation operators, /- and f, , 
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which satisfy the usual fermionic anticommutation rules, 

{fi ,/,•} = {fl /}} = 0, {fi , /]} = Sij. This leads to a 
model that describes free spinless fermions, 



h.c. 



wY,^-L/2fn tl (5) 



where n, = fjf i represents the number of fermions on 
site i. Because the model (J5J) is a quadratic form of j\ and 
/j , it can be solved by a simple numerical diagonalization 
of the L x L matrix. Denoting by with k G [1,L] the 
eigenvalues of this matrix, the partition function ([2]) takes 
the form 



n 
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The entropy is defined as S — — fcsTr V\nT> with the 



density matrix T> 



1 



Working in a sys- 



tem of units where the Boltzmann constant kg = 1 
and using the properties of the density matrix, it fol- 
lows that S = InZ + - fi^(N). Substituting 
d 

(H) — /i(A/") = — — InZ and using expression for 
Z, the entropy takes the form 



S(J3,n) 



L 

E 

fc=i 



In 1 



e /3(e fc -/i) _|_ l 



(7) 



The average number of particles N is obtained by sum- 



ine average nuniDcr or particles i\ 
ming the Fcrmi-Dirac distribution, 



iV(^M) = E' 



1 
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Fig. Q] shows the entropy as a function of temper- 
ature for both DC and ODC cases. The chemical po- 
tential fi is adjusted such that the average number of 
particles (jSJ) remains constant (N = 50). An interest- 
ing feature is that the two curves cross at a temperature 
T c , and that below T c the entropy of the ODC system is 
greater than the entropy of the DC system. Thus, if the 
initial temperature Td c is below T c , then the final tem- 
perature T D dc is lower when switching adiabatically from 
DC to ODC. 

Next we generalize our discussion by calculating, for 
a fixed number of particles N, the critical temperature 
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FIG. 1. (Color online) The entropy as a function of temper- 
ature for 50 hard-core bosons, in the DC case (circles) and 
in the ODC case (triangles). There exists a critical tempera- 
ture T c where the two curves cross. If the initial temperature 
Td c is below T c , then the conservation of the entropy when 
switching adiabatically from DC to ODC implies that the final 
temperature T d c is lower. 

T c below which the ODC method produces a tempera- 
ture lower than the DC method when the confinement is 
switch adiabatically. In order to determine T c for given 
parameters tf? and W for the conventional DC system, 
and i° dc for the ODC model, one needs to solve for each 
value of /V a system of three coupled non-linear equa- 
tions, 

SodcWc Hodc) — SddPc-, Mdc) 

N odc (/3 C) node) = N (9) 
N dc (p c , fide) = N 

where S odc (S dc ) is given by Eq. © and N odc (N dc ) is 
given by Eq. ©, with fj, = fi odc (p = fj, dc ), and P = C . 
The first equation corresponds to the conservation of the 
entropy when switching from DC to ODC, and the two 
others correspond to the conservation of the number of 
particles. Solving this system of equations determines the 
critical inverse temperature /3 C — 1/T C , and the chemical 
potentials \i odc and \i dc that give the desired number of 
particles N. 

For this purpose, we define an error function: 

£(/3 c , Mode, Mdc) = (Sodc(Pc, Hodc) ~ S dc (/3 C , /X<jc)) 

+ {N odc (f3 c ,ti odc ) -N) 2 

+ (N dc {f3 c ,ti dc )-N) 2 . (10) 

By construction, the solution of Eq. (JSJ minimizes this 
error function. Starting with an initial guess for j3 c , [i odc , 
and [i dc , we calculate the error £(/3 c , Mode, Mdc) and its 
gradient V£ = (d£/df3 c ,d£/dn odc ,d£/dii dc ). Writing 
the initial guess as a vector, f — (/3 C , AWc, Mdc), we per- 
form a correction Ar by following the opposite direction 
of the gradient V£. Then we iterate until convergence. 



FIG. 2. (Color online) The critical temperature T c and the 
DC isotherms as functions of the number of particles TV. At 
T — T c , there is no change in temperature when switching 
adiabatically from DC to ODC. Below T c (blue region), the 
ODC method gives a temperature T d c that is lower than the 
temperature Tdc obtained with the DC method. Above T c 
(yellow region), it is the DC method that gives the lowest 
temperature. For example, the point P corresponds to a sys- 
tem with 34 particles, an initial temperature Tdc = 3, and a 
final temperature T d c ~ 1-5. 

Fig. [5] shows the critical temperature T c and the DC 
isotherms as functions of N. For a given number of parti- 
cles and an initial temperature T = T C (N), the ODC and 
DC systems have the same temperature T odc — T dc when 
the confinement is switch adiabatically. Below (above) 
T c , the ODC system has a temperature T odc that is lower 
(higher) than the temperature T dc of the DC system. The 
point P illustrates how the figure should be read: For a 
system with 34 particles and an initial DC temperature 
T dc — 3, the final ODC temperature is T odc w 1.5. Note 
that T c vanishes when N=L=70. The resulting Mott 
phase found in the ODC case always has lower entropy 
than the mixed phases found in the DC case. This will 
be discussed in greater detail in the next section. 



IV. QUANTUM MONTE CARLO ALGORITHM 
AND THE ENTROPY 

For the treatment of soft-core interactions, we perform 
QMC simulations using the Stochastic Green Function 
(SGF) algorithm [11] with tunable directionality [Hj]. Al- 
though this algorithm was developed for the canonical en- 
semble, a trivial extension [Tfj| allows us to simulate the 
grand-canonical ensemble. We propose a new method to 
measure the entropy by taking advantage of the grand- 
canonical ensemble. Our thermodynamic control param- 
eters are the temperature T, the volume V (number of 
sites L), and the chemical potential fi. Unlike the analyt- 
ical hard-core case, a direct measurement of the entropy 
is not possible with a single QMC simulation because the 
value of Z is unknown. However it is still possible to eval- 
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FIG. 3. (Color online) The entropy S and the thermal suscep- 
tibility x m the hard-core case. Comparison between results 
obtained with exact diagonalization using Eq. © and QMC 
results with the SGF algorithm and Eq. ()15|) . Two different 
temperatures are considered, T = 3 and T = 0.25, for both 
DC and ODC systems. 
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FIG. 4. (Color online) The entropy as a function of tem- 
perature for 50 soft-core bosons with U — 8, in the DC case 
(circles) and in the ODC case (triangles). As in the hard- 
core case (Fig[T|, there exists a critical temperature T c below 
which the entropy of the ODC system is higher than the en- 
tropy of the DC system, thus making the ODC method more 
efficient than the DC method for producing low temperatures. 



uate the entropy with a set of QMC simulations. For this 
purpose, we define the thermal susceptibility Xth by the 
response of the number of particles N to an infinitesimal 
change of the temperature T: 

By substituting N = -§Tr J\f e ^ p( ^^ in expres- 
sion (jlll) . we get an expression for the thermal suscepti- 
bility that can be directly measured in our simulations: 



xth = P- 



(12) 



is defined as P — ——— 
oV 



Considering the energy E = {%.) and the associated dif- 
ferential dE = TdS — PdV + fidN, where the pressure P 

, and performing a Legendre 

S,N 

transformation over the variables S and N, we can define 
the grand-canonical potential Q that depends only on our 
natural variables, Q(T,V,fi) = E — TS — fiN = -PV. 
Its differential takes the form 

dQ = -SdT - PdV - Ndfi. (13) 

We can then extract a useful Maxwell relation, 



dS_ 
<9/i 



V,T 



dN 



(14) 



so the entropy can be easily obtained by integrating the 
thermal susceptibility over the chemical potential and 
keeping the temperature and the volume constant, 



(15) 




FIG. 5. (Color online) The entropy as a function of the 
inverse onsite interaction 1/U for 50 soft-core bosons at T = 1, 
in the DC (circles) and ODC cases (triangles). The entropy of 
the ODC system remains greater than the entropy of the DC 
system for any value of the onsite repulsion U, showing that 
our suggested cooling method works deep into the soft-core 
case. 



where /io is the critical value of the chemical potential 
below which the average number of particles iV and the 
thermal susceptibility xth are vanishing. 

In order to check the reliability of Eq. (fT5|) . we show 
on Fig. [3] a comparison of the entropy of the hard-core 
case obtained with the SGF algorithm by integrating the 
thermal susceptibility (|12[) . and the entropy computed 
with Eq. 0. The agreement is good for both DC and 
ODC systems at high (T = 3) and low temperatures 
(T = 0.25). 

We now release the hard-core constraint and set the 
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FIG. 6. (Color online) The entropy as a function of the 
inverse onsite interaction 1/U for 70 soft-core bosons at T — 1, 
in the DC case (circles) and in the ODC case (triangles). For 
the ODC case the filling is commensurate and the system 
forms a pure Mott phase that cannot carry any entropy. Thus 
ODC cannot be used as a cooling method for commensurate 
fillings. 



onsite repulsion [7 = 8. Fig. @] shows the entropies for 
the DC and the ODC models as functions of tempera- 
ture for N = 50. The curves differ from the hard-core 
case only quantitatively, not qualitatively, showing that 
the method of cooling by switching from DC to ODC 
still works. Moreover, one notices that the critical tem- 
perature T c 3.5 is higher than in the hard-core case 
(T c sa 2.65) which makes easier to access the regime 
where ODC is more efficient than DC. 

Further, we extend our soft-core results to different 
values of the onsite repulsion. Fig. [5] shows the entropy 
for the DC and the ODC models as function of the inverse 
onsite repulsion 1/U for N = 50 and T — 1.0. The curves 
show that the entropy of the ODC model is above the one 
of the DC model for any value of U. Thus, for this filling, 
the ODC method produce temperatures lower than the 
DC method for any value of the onsite repulsion. When 
U is large, the results match with those obtained in the 
preceding section for the hard-core case (Fig. [1]). 

However, the situation is different for N = 70 as Fig. [6] 
illustrates. At this integer filling, the entropy of the ODC 
model, which vanishes in the large U limit, intersects 
the curve for the DC model. In this regime, the ODC 
model exhibits a pure Mott phase, hence with zero en- 
tropy. However the phase of the DC model has Mott 
regions coexisting with superfluid regions, so the entropy 
remains finite. Thus, ODC cannot be used to cool the 



system in this region. 

Concerning the experimental realization of our model, 
a holographic technique recently developed [TtJ can be 
used to build the optical lattice with off-diagonal con- 
finement. Using this method, an off-diagonal trap can be 
superposed to an existing diagonal trap. Then the diag- 
onal trap can be turned off. The switching between the 
two traps can be in principle very fast, however the tech- 
nical details of how this will work go beyond the scope of 
the present manuscript and must be developed by exper- 
imentalists. Nevertheless, a qualitative analysis reveals 
that three time scales must be considered. The time 
scale r m of the model system or roughly r m = 1/t (in 
units where h = 1), the time scale of the experiment 
r e , and the time scale r c which describes the coupling 
of the model system to its environment which includes 
the effects of the laser heating, evaporation, etc. In our 
proposal, it is important that the trap is adiabatically 
switch on the experimental time scale, but not on the 
time scale which describes the coupling of the trap to its 
environment, so that r m << r e << r c . 

V. CONCLUSION 

In this manuscript we propose that the adiabatic 
switch from the DC to the ODC method can produce 
lower temperatures for a wide range of initial tempera- 
tures and system parameters. In the hard-core limit, we 
determine the critical temperature T c for which the two 
methods have the same entropy. Below (above) T c and at 
constant entropy, the ODC method leads to temperatures 
that are lower (higher) than with the DC method. In or- 
der to extend our results to the soft-core case, we propose 
a simple method for evaluating the entropy with QMC, 
by measuring the thermal susceptibility xth i n the grand- 
canonical ensemble and integrating it over the chemical 
potential /i. Then we make use of the SGF algorithm [l4| 
with tunable directionality [lj|, and show that the soft- 
core results are qualitatively the same as in the hard-core 
case. 
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